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Abstract An experiment in which the third Stokes parameter thermal emission from a 
periodic water surface was measured is documented. This parameter is shown to be related 
to the direction of periodicity of the periodic surface and to approach brightnesses of up 
to 30 K at X band for the surface used in the experiment. The surface actually analyzed 
was a ‘two-layer’ periodic surface; the theory of thermal emission from such a surface is 
derived and the theoretical results are found to be in good agreement with the experimental 
measurements. These results further the idea of using the third Stokes parameter emission 
as an indicator of wind direction over the ocean. 

1 Introduction 


Recent theoretical works have suggested the potential of passive polarimetry in the remote 
sensing of geophysical media [1-3]. These works indicate that the third Stokes parameter, 

Ub , of the thermal emission may become large for azimuthally asymmetric fields of ob- 
servation. In [3], values of Ub as high as 40 K were measured from a periodic triangular 
soil surface at 10 GHz. In order to demonstrate the existence of the Ub parameter for 
a water surface, and thus to investigate the potential applicability of passive polarimetry 
to ocean remote sensing, observations of a periodic water surface were performed at the 
Cold Regions Research and Engineering Laboratory (CRREL), Hanover, NH in Septem- 
ber of 1992. This paper documents those observations and the theory developed for the 

“two-layer” periodic surface actually analyzed in the experiment. 
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2 Theory of Polarimetry 


In passive polarimetry, brightness temperatures corresponding to all four modified Stokes 
parameters are measured. The brightness temperature Stokes vector is defined as 



\Ih 


(EhEZ) 

1 

u 

1 

(E.EZ) 

c 

u 

~vC 

2Rc{E,E' h ) 


. V . 


2lm(E v E^) 


( 1 ) 


In the above equation, Eh and E v are the emitted electric fields received from the horizontal 
an d vertical polarization channels of the radiometer, 77 is the characteristic impedance, 
and C = K/ A 2 with K denoting Boltzmann’s constant, A the wavelength. The first two 
parameters of the brightness temperature Stokes vector correspond to the received powers 
for horizontal and vertical polarizations, respectively. The third and fourth parameters 
correspond to the complex correlation between the electric fields received by the horizontal 
and vertical channels. We label the four parameters Tjsh, Tbv> Ub> and Vb respectively in 
this paper. 

It is shown in [2] that the third and fourth Stokes parameters may be related to the 
brightness temperature in a 45 degree linearly polarized measurement ( Tb p ) and a right- 
hand-circularly polarized measurement ( Tb t ) as follows: 


U B = 2 Tbp-Tbh-Tbv (2) 

Vb — 2Tbt — Tbh — Tb v (3) 


Thus, to measure all four parameters of the Stokes vector, the brightness temperatures in 
horizontal, vertical, 45 degree linear, and right-hand-circular polarizations are first mea- 
sured, and the above equations are used to obtain Ub and Vb- 

In the passive remote sensing of rough surfaces, the parameter that is actually of 
interest is the emissivity, which relates the brightness temperature emitted by an object 
to its actual physical temperature, under the assumption that the object is at a constant 
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physical temperature and that the emission from the object is the only source of brightness: 

Tbo. = e a (Q, <t>)Tph y , (4) 

In the above equation, the subscript a refers to the polarization of the brightness tempera- 
ture, 0 to the polar observation angle, and (f> to the azimuthal observation angle. Through 
the principles of energy conservation and reciprocity, Kirchhoff ’s Law relates this emissivity 
to the reflectivity of the surface [4]: 

ea(0,<f>) = 1 ~r a {0><t>) (5) 

The reflectivity r a (0 , <f>) for the given incident polarization a is defined as the fraction 
of the power incident from direction ( 9 , <fi) that is rescattered and can be evaluated by 
integrating the bistatic scattering coefficients 'yi )a (0,<f>',0 l ,<f> > ) over all scattering angles in 
the upper hemisphere and summing the results of both orthogonal scattering polarizations. 

r a (9, <f>) = £ J*' 2 d9' sin O' j*' d<f>' lba {9 , 0; 9', <f>') (6) 

b 

In the expression of the bistatic scattering coefficient, (0 t <f>) and (9',<f)') represent the 
incident and the scattered directions, respectively, and the subscripts a and b represent 
the polarizations of the incident and the scattered waves, respectively. 

Thus, to calculate the fully polarimetric emission vector, the bistatic scattering co- 
efficient for each of four polarizations is first calculated and integrated over the upper 
hemisphere to obtain the reflectivity for that particular polarization. Multiplication of the 
corresponding emissivity by the physical temperature of the object under view yields the 
brightness temperature for this polarization. The fully polarimetric brightness vector is 
then calculated as described previously. 

3 Two-Layer Periodic Surface Theory 

In order to form a periodic water surface, a thin sinusoidal sheet of fiberglass was placed 
on top of a flat water surface. After air bubbles trapped underneath were removed, a 
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“two-layer” periodic surface resulted. To examine the effects of this layer on the measured 
brightness temperatures, a theoretical model for thermal emission from a two-layer periodic 
surface was derived and implemented. 

The problem of calculating thermal emission from a “single-layer” periodic surface has 
been solved by several researchers using many different techniques. Two of the methods 
used in the single-layer case are extended for the two-layer case in this analysis. The results 
from these two ‘exact’ methods are then compared to determine which would apply best 
to the cases measured in the experiment. For both of these methods, the calculations are 
performed to determine the total power reflected in the reciprocal active problem; this 
total reflectivity is then used to obtain the emissivity as described previously. 

The extended boundary condition (EBC) method described by Chuang and Kong [5] 
for calculating scattering from a single-layer periodic surface is the first method to be 
applied. In this method, Huygens’ principle is applied at the surface so that the scattered 
field can be obtained once the surface fields are known. The Huygens’ integral equation is 
solved by expanding the unknown surface fields into Fourier series and forming a truncated 
impedance matrix which is inverted to obtain the surface fields. This method is known as 
the extended boundary condition method because the calculations enforce the requirement 
that the Huygens’ integral which equals the field value within the region of interest must be 
zero outside that region. Although the extended boundary condition method is exact, the 
numerical formulation becomes ill conditioned if the surface height-to-period ratio becomes 
too large [10]. 

To overcome this problem, a second method, the method of moments [9], can be used 
to perform the calculations. In this formulation, the same Huygens’ principle equation is 
used, but the requirement that the Huygens’ integral in the region not of interest be zero 
is not enforced. Instead, the integral equation is “tested” at a discrete number of points 
along each surface where the equation is forced to hold. As the number of testing points 
is increased, the calculated scattered fields converge to the exact results. The method of 
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moments can be used in cases where the extended boundary condition method fails. It 
requires, however, more computation time than the EBC. 

3.1 General Formulation 

Consider a plane wave incident upon a two-layer periodic surface described by fi(x) = 
fi(x + P), U (x) = f 2 (x + P ) with fi representing the upper surface, f 2 the lower surface 
and P denoting the period of the surface in the x direction (see figure 1). The electric field 
of the incident wave is given by 

Ei = iiE 0 e ilcir (7) 

where ki denotes the incident wave vector and is equal to xk x i + ykyi — zk z i and e; is the 
polarization of the electric field vector. 

Since the structure is uniform in the y direction, all the field components in regions 0,1, 
and 2 will be phase matched to the same exp (ikyiy) dependence. Thus, we can replace all 
the terms in Maxwell’s equations with ikyi as is commonly done in waveguide theory 
(usually with the z dependence by convention.) Subsequent equations will thus have the 
exp (ik^y) dependence removed. 

Maxwell’s equations can now be simplified so that the x and z components of the 
electric and magnetic fields can be expressed as functions of the y components of the fields 
as follows: 

Ej *00 = , 2 \ 2 l k vi V *- E iv( f ) + x Hjy(f)] (8) 

Kj Kyi 

Sj $ (T) = — j 2 ~ X Ejy^f)] (9) 

Kj Kyi 

where j = 0,1,2 signifies regions 0,1, or 2, respectively, V, is the gradient operator that is 
transverse to the y direction 

v -=4 + 4 (io) 

and a = ( x,z ) indicates the transverse components of the fields for region j. The y 


5 


components of the electric and magnetic fields now satisfy the partial differential equations 


(V, ! + kj* - = 0 


(ii) 


where Aj y = Ej v , Hj y . Since V, operates only in a two-dimensional space, the y components 
of the electric and magnetic fields satisfy a two-dimensional wave equation. The Green’s 
function for such an equation is given by 


<?,(/>., p'.) = - £|) (12) 

= 4- [ dk x y— x exp [ik x (x — x') + ikj x \z — z'|] (13) 

47t y-oo kj z 

where j = 0,1,2, 

kj, = (k/ - (14) 

ki, = M - kl ( 15 ) 

and 

p, — xx + zz (16) 


Applying the scalar Green’s function form of Huygens’ principle to the top periodic 
surface and region 0 gives 
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where the y dependence has been suppressed and 
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The above surface integral is over an infinite domain. However, the periodic properties 
of the surface fields can be used to reduce the integration domain to a single period. The 
surface fields have the property 


oi(p # -f xtiP ) = u>(p J )e ,fcc * ni> 


(19) 
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where u(p t ) can be Ej y (p s ), Hj y (p,) or their normal derivatives, n is an integer, and P is 
the period of the surface. 

Following Chuang and Kong [5], a periodic Green’s function can be derived using this 
property. This Green’s function is 

Gjp{p„p t ) = 0 ., 1 p X) ex P [ikj.a jn (x ~ x') + ik t /3j n \z - z'\) (20) 

llKj t r n Pj n 
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Now using the periodic Green’s function in the scalar Huygens’ principle formulation 
for each of the three regions yields 
Region zero 
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where (5"/) indicates that the integration is performed over surface l. 

Once the surface fields are determined, the fields in regions 0,1, and 2 can be obtained 
using these relations. At present, there are sixteen unknown vectors to be obtained: the y 
components of the electric and magnetic fields at surface one for regions zero and one, the 
y components of the electric and magnetic fields at surface two for regions one and two, 
and the normal derivatives of the y component of the electric and magnetic fields for each 
of these surface-region combinations. Electric and magnetic field boundary conditions can 
be used at both surfaces to reduce the number of unknown vectors to eight. 


3.2 Boundary Conditions 


At both surfaces, the y component of the electric and magnetic fields are tangential, and 
therefore must be continuous: 
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Eoy = Ely z = fi(x) 

(29) 

$ 

II 

(30) 

E 2y = E 2y z = f 2 (x) 

(31) 

n 2y = E 2y 

(32) 

The components of the fields along the surface profile are 
continuous. These relations can be written as: 

also tangential and must be 

ni x E 0 , = n 1 x E u z = /i(x) 

(33) 

ni x Ho. = fix X Hu 

(34) 

n 2 x E u = n 2 x E 2 , z = f 2 (x) 

(35) 

fi 2 x Hi. = fi 2 x H 2 . 

(36) 


Note that both the normal and the s component fields lie in the x — z plane, so that 
hj x A, is a y directed vector whose amplitude corresponds to the tangential component 
of the field along the surface profile. Following Chuang and Kong [5], the second set of 
boundary conditions can be written as: 

y{ni • V,2?o„) = c 0 hi x + c 2 y(h x • V.E 2y ) z = fi(x) (37) 

y(n 1 • V,H 0y ) = -d 0 ni x V t E ly + d 2 y(fii ■ V.H ly ) (38) 


where 
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y(h 2 • V.H ly ) = -d x n 2 x V ,E 2y + d z y{h 2 • V.H 2y ) 
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fio k 2 u /zx kl 

These eight boundary conditions reduce the sixteen unknown vectors to eight. The 
remaining eight unknowns are then solved by expanding the surface fields in the Huygens’ 
integral equations into a sum of unknown amplitude basis functions. The functions used 
in this expansion and the equations solved differ for the EBC and Moment methods; each 
method will be treated separately beginning with the EBC. 

3.3 Extended Boundary Condition Method 

H we substitute the expression for the periodic Green’s function into the Huygens’ Principle 
integral equations and factor out the unprimed exponential dependence, the equations 
reduce to: 
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where the |z - z'| expression in the periodic Green’s function has been replaced by z - z' 
when the observation point is above the surface, z 1 — z when the observation point is below 
the surface. In the above equations, 
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and and c^\ have the same form as their electric field counterparts with Ej y replaced 
by Ejy. We now recognize that the scattered field in region 0 is composed of a discrete set 
of plane waves traveling in the directions indicated by These plane waves are known 
as the Floquet modes; the amplitude of the nth reflected mode is determined by 6 0 i n * 

For the extended boundary condition method, we use the equations which involve the 
region of no interest for each region. These are equations (46), (48), (49), (51), (52), (54), (56), 
and (58). These eight equations will uniquely determine the eight unknown surface field 
vectors. Once the surface fields are determined, the remaining equations are used to 
calculate the scattered field in each of the three regions. 

Following Chuang and Kong [5], a Fourier series expansion is used for each of the 
unknown field quantities as follows: 
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Substituting the above expressions into the eight equations and using the property 


dan, x V t H jy (x,z{x)) - ydx ^ + dx dz 

leads to the coupled matrix equation of Appendix A. Truncating the Fourier series of the 
surface field unknowns results in a finite matrix equation which can be solved to obtain 
the surface field unknowns. The upward going field amplitudes are then calculated to be 
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and the downward going field amplitudes in region 2 are 


(71) 

(72) 


U22 = E 2 (3 B +D 2 a B (73) 

aff =I 2 “6 b +^ 2 V (74) 


The matrices in the above equations are as described in Appendix A. The above coefficients 
are the amplitudes of the y components of the scattered and transmitted electric and 
magnetic fields. The x and z components can be calculated from equations (8) and (9). 
The scattered plane waves in region zero propagate along the directions determined by 
kQ n , while the transmitted plane waves in region two travel along the directions of k 2n . 

Once the amplitudes of the scattered modes are calculated, the emissivity of the surface 
can be obtained by integrating the toted power reflected over the upper hemisphere and 
then applying Kirchhoff’s law to obtain the emissivity corresponding to the polarization 
of the incident field. For the n th propagating scattered Floquet mode, the reflectivity is 
given by power transmitted in the z direction divided by the incident power and is found 
in [5] to be: 


Ifeoml 2 + Mff n | 2 

laoiol 2 + l^oaffol 2 
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where the a 0 io and ao^ are as defined in Appendix A. The brightness temperature of the 
surface is then given by 

Tsa = ( 1 — y> n ) (76) 

n 

where the sum is over the propagating modes reflected for an incident wave polarized in 
the a direction. 


3.4 Method of Moments 


A second method for solving the Huygens’ integral equations for a two-layer periodic surface 
is the method of moments. The use of the method of moments for computing scattering 
from a single-layer periodic surface has been presented by Veysoglu et al [2]; these methods 
will be extended for the two-layer periodic surface in this section. 

We begin with the integral equations presented in the general formulation section and 
the periodic Green’s function. We will also use a second form of the periodic Green’s 
function in the calculations 

Gjp(p,p') = J X - (x' -f rap )} 2 + (z- z') 2 ) exp (ik 0 x mp) (77) 

^ m= — OO 

in addition to the sum of exponentials form used in the EBC method. The eight surface 
field unknowns in the integral equations are now expanded into sums of unknown amplitude 
“pulse” functions as follows: 
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where 
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and N is the total number of basis (pulse) functions used. 

Instead of solving the equations with right-hand-sides of zero, as was done in the EBC 
method, the method of moments uses the equations with right-hand-sides corresponding 
to the fields in the region of interest. These six equations are “tested” at a discrete set of 
points along the surface upon which the integral is performed where the integral equations 
are forced to hold. The missing two equations needed to match the eight unknown functions 
are obtained from testing the integral equations for region one on both surfaces one an d 
two. The testing points are chosen to lie in the center of the pulse basis functions described 
above. The integral equations now reduce to the matrix equation given in the Appendix 
B. 

Note that this method will require the integration of the singularity in the periodic 
Green’s function when the testing point and integration range overlap. For these “self” 
terms of the matrices, an asymptotic expression for the Hankel function of a small argument 
is used and analytically integrated. The integral expression proposed by Veysoglu et al [2] 
is also used to speed up the convergence of the evaluation of the Hankel function sum in 
the periodic Green’s function. After solving the matrix equation to determine the surface 
fields, the coefficients of the reflected and transmitted Floquet modes can be calculated. 
If the reflected fields in region zero are written as 
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and the transmitted fields in region two as 
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the coefficients are given in terms of the surface fields by 


b = -i_ f 

n 2ik 0 .P Jr 


o.P JP{S 1) 


dS 


e-'Kn-P 


fil • V.Eoy(p) 


~Eo y {p)fll • V # 




bW = - 1 - [ 
n 2ik 0t P Jr 


P(51) 


dS 


e-'hn-P _ 


Pn 


n x • V t H 0y (p) 


-H 0y (p)h 1 • V, 




(91) 


(92) 


A n 


l —r 

2 ik 2 ,P JP(S2 ) 


dS 


’-'hnP 


-h 2 • V,E 2y {p) 


-E 2y {p)n 2 • V,- 




& 




= -J-/ 

2ik 2 ,P Jp{S2 ) 


dS 




/?« 


n 2 • V t H 2y (p) 


-H 2y {p)h 2 • V, 




(93) 


(94) 


where is defined as in the EBC section. 

As with the EBC, once the mode amplitudes are determined, the emissivity can be 
found by integrating the total power reflected over the upper hemisphere and then applying 
Kirchhoff ’s Law. Note that the method of moments does not suffer from the problems of the 
EBC when the height- to-period ratio becomes too large. A larger number of basis functions 
is required and the calculations take longer to perform, but the method ultimately yields 
an accurate result. 
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3.5 Comparison of Methods 

Both the EBC and MOM methods were used to perform a theoretical analysis of the sur- 
face in the experiment. Figure 2 is a plot of the four polarimetric brightness temperatures 
calculated by both the EBC and MOM. Also shown in these plots are the corresponding 
brightness temperatures for a periodic water surface without the fiberglass layer on top. 
These calculations are for the experimental parameters described in the next section as- 
suming a physical temperature of 300 K. The EBC and MOM results axe seen to agree, 
indicating that the surface slope limitations of the EBC axe not being exceeded and that 
the numerical resolution of the method of moments code is adequate. The more efficient 
EBC method was used for the rest of the theoretical calculations presented. It is also 
seen that the effect of the layer is to increase Tbh> Tj 3 „, and Tb p by approximately 3 OAT 
over the single-layer brightnesses, but the effect on Ub is small due to cancellation of this 
increase when the brightnesses in the three polarizations are used to compute Ub . Finally, 
the characteristic “jumps” in brightness temperature for an infinite periodic surface due 
to the transition of a Floquet mode from propagating to non-propagating are seen in the 
Tbk curve. Note that the exact shape of the curve may not be the same as that of Figure 
2, as the curve was linearly interpolated between the calculated points every 15 degrees in 
azimuth. 

4 Experiment Setup 

To form the periodic water surface, a sinusoidal sheet of fiberglass 1.32 m x 1.22 m x 
0.002 m thick was placed on top of a 1.83 m side square pool of fresh water of 0.23 m 
depth. The height and period of the sinusoidal corrugation were 1.4 cm peak-to-peak and 
6.8 cm respectively. To prevent water from flowing over the edges, a styrofoam border 
approximately 3 cm high was glued around the surface. 

The temperature of the water was measured around the surface in the experiment and 
found to be 22.5° ± 1.5° C throughout the experiment. As fresh water was used, this 
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corresponds to an average dielectric constant of 62 + 132 at 10 GHz over this temperature 
range [7]. The dielectric constant of the fiberglass layer was measured using a network 
analyzer technique [6] from 12 to 18 GHz and found to be 3.1 ± 0.1 + 1(0.05 ± 0.05) over 
this entire band. It was assumed that this dielectric constant was the same at 10 GHz. 

The radiometer used in the experiment is the same as that used in [3]. It operates 
at 10 GHz with a 200-MHz RF bandwidth, —3 mV/K sensitivity, ±1 K precision, and 
a pyramidal (standard gain) horn of 30° beamwidth. The radiometer was mounted on a 
tripod at an elevation of 1.7 m above the pool surface and directed toward the surface 
along a direction determined by the azimuthal angle, <f> } and the polar angle, 6 (see Figure 
1). The radiometer was calibrated by viewing an absorber whose physical temperature 
Wets known for the hot load, and by viewing the sky at a set of angles for the cold load. 
Since the measured sky voltage can be expressed as a function of polar angle, the points 
measured can be extrapolated to obtain a radiometer voltage corresponding to brightness 
temperature 3 K. This voltage is then used as the cold reference for the calibration. 

Observations of the periodic surface were made as a function of azimuthal angle, which 
was varied by rotating the periodic surface in the pool so that the background noise contri- 
bution from the antenna sidelobes would remain constant through a set of measurements. 
The horn of this single polarization radiometer was rotated to produce measurements of 
Tbaj Tbv, and Tb p . Measurements of Tb t were not possible with this radiometer; however, 
the Vb brightness was not predicted to be significant. 

5 Estimation of Noise Contribution 

The contribution of background noise to the measured antenna temperature for the periodic 
surface can be large: noise from the antenna sidelobes and from reflections of sky bright- 
ness off the surface both contribute to error in the measurements. Since this background 
noise could not be easily eliminated from the experiment, estimates of the background 
contribution were made so that the noise could later be removed. The contribution of the 
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background noise was estimated by measuring the brightness temperature of a flat water 
surface alone and by measuring the brightness temperature of a reflector in the pool which 
was the same size as the periodic surface. The theoretical results for both of these cases 
are known; thus, these measurements form a sort of two-point calibration for the periodic 
surface when it is in the pool. 

Consider the case of the observation of the water in the pool alone. The measured 
brightness temperature, TfiatM > consists of contributions from the water surface, T// at , 
from the ground beside the pool and the sides of the pool, T,i<Uy smd from reflections of 
the sky temperature, T t k y , off the pool. The contribution of the reflections of the sky 
temperature from the ground beside the pool is neglected. If we assume that there is no 
variation in these temperatures with the incidence angle, then the relationship between 
these variables can be written as: 

TfiatM = Tflat( 1 — ft) + T,ky{ 1 — ft){ 1 — Tflat/Tp wa t ) + T,ideft (95) 

where f, represents the fraction of the antenna pattern that does not view the pool and 
Tp wa t is the physical temperature of the water in the pool. Measured values of f, ranged 
from 3 to 18 percent through the experiment. Both T s k y and T t ide were obtained from the 
experiment at angles corresponding to the specular reflection for the sky temperature and 
direct incidence for the ground temperature. 

Next consider the case of a reflector in the pool the same size and at the same loca- 
tion as the periodic surface. The brightness temperature measured, T re fiM , consists of the 
reflection of the sky temperature over the reflector surface, the flat water surface temper- 
ature over the remainder of the pool and the reflection of the sky temperature from this 
portion of the surface, and the ground contribution. This can be expressed as: 

TrefiM — T,k y fr + T, jfc y (l — f T ~ /«)( 1 — Tfl a t/Tp wa t) 

+Tflat(l — f T — ft) + T tide ft 
= TfUtM + Tflatf r (Ttk y /Tp W at ~ 1) 
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(96) 

(97) 


where f T is the fraction of the antenna pattern that views the reflector surface, which was 
found to range from 50 to 74 percent through the experiment. A further assumption that 
the physical temperature of the water and the background sky temperature are the same 
in the flat surface and reflector measurements is made in the second equality. 

Finally, consider the case of the periodic surface in the pool. The measured temperature, T, ut jmi 
consists of contribution from the periodic surface, the flat water in the pool, sky reflections 
off the pool area, and the side contribution. The reflections of sky temperature off the pe- 
riodic surface are determined by its Floquet modes; however, since the sky temperature 
is assumed uniform, this reduces to the same case as that of a specular reflector. The 
measured temperature can be expressed as: 

TturfM = Tnrffr + T tk J r (l - T turf /T Pwat ) + T //at (l - - /,) 

+T jfc „(l — Tflat/Tp wat )( 1 — fr — ft) + Ttideft (98) 

= T re flM + T tur ff r ( 1 - T, k yfTp wat ) (99) 

These three equations can be used to solve for T, ur /, the periodic surface brightness 
temperature averaged over the fraction of the antenna pattern on the surface: 

t.„, = _ T % — (r. WM - (ioo) 

■LflatM ~ 1 refl M 

This brightness temperature has the effects of both the reflections of the sky temperature 
off the surface and the antenna field of view outside of the periodic surface removed. 

This calibration procedure increases the effect of the radiometer measurement uncer- 
tainty in the final calculated T, ur f to approximately 3 K. The assumptions made in the 
above equations also contribute to the error in the experiment. A theoretical study in 
which the brightness contributions of the water around the surface and the sky tempera- 
ture were estimated by integrating the angularly varying T/; at and T, jt y over a theoretical 
antenna pattern indicated that the use of the specular Tji at and T t k y values were valid to 
within 2 K. A bias in both of these values would produce a systematic error of less than 
5 K. One remaining possible source of error is the possibility of the close-to-the-horizon 
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Floquet modes of the periodic surface imaging something with a higher brightness temper- 
ature than the sky. Note that the Ub parameter is fairly insensitive to any systematic error 
in the calibration used. Since systematic error will tend to affect all three polarizations in 
a similar manner, the effects tend to cancel out, leaving Ub unaffected. 

The results presented in Section 6 for the measurements of periodic surface brightness 
temperatures were obtained using the above technique. 

6 Experimental Results 

Table 1 presents the measurement results at Q — 20°. Shown are the measured emissivities 
eA,e v , and e p . The emissivity of the third Stokes brightness parameter is then obtained as 
eu = — 2e p . This formula is opposite in sign to that previously given due to the 

fact that the radiometer polarization basis was actually rotated to —45° in the experiment. 
It is observed from the experiment that, in general, the cu value approaches a maximum 
absolute value at <j> = —45° and approaches zero at <f> = —90° and <f> = 0°. This shows 
that the Ub parameter is sensitive to the azimuthal direction of a periodic water surface. 
These trends are also seen in Table 2, which presents the measurement results at 6 = 30°. 

7 Comparison with Theory 

Figures 2 and 3 are plots of the measurement results and the theoretically predicted bright- 
ness temperatures for polar angles 20° and 30°, respectively, normalized to a constant 
physical temperature of 22.5 C. The theoretical values plotted were obtained using the 
two-layer version of the Extended Boundary Condition method. These theoretical results 
were also averaged noncoherently over 81 points in a three-dimensional calculated antenna 
pattern for the pyramidal horn used in the experiment, which tends to smooth out the 
abrupt changes obtained for plane wave incidence in section 3. This pattern was obtained 
from [8], using the measured antenna dimensions of 5.7 cm by 7 cm for the aperture di- 
mension and 10 cm and 8.9 cm for the corresponding flare lengths of the horn. Good 
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general agreement is observed between the theory and experiment, indicating that the 
background noise removed procedure was reasonable. However, there is a bias of the ex- 
perimental results above the theoretical predictions of approximately 10 K for Tb v 
and Tb p and 6 K for Ub- One possible source of this error is the permittivity value for the 
fiberglass surface used in the experiment, which could be slightly higher at 10 GHz than 
the measured value at 12.4 GHz. A increase in this permittivity to 3.2 -f- zO.l increases the 
theoretical predictions by approximately 5K or 0.017 emissivity. A second possible source 
of the error is the use of the specular values for T/i at and T,k y as discussed in section 5. 
This could also contribute approximately 5 K or 0.017 emissivity variation. Finally, the 
imaging of a brightness source other than the sky by the Floquet modes of the periodic 
surface would tend to increase the experimental brightness temperatures above the theo- 
retical predictions. However, the area in which the experiment was performed was clear 
to about ten degrees above the horizon. Only the azimuthal angles —90, —15, 0, and 15 
are theoretically predicted to have modes this close to the horizon. 

8 Conclusion 

We have demonstrated in this paper that the Ub parameter exists for a two-layer periodic 
water surface and can approach brightnesses exceeding 30 K at X band. We have also 
shown that the Ub parameter is fairly insensitive to the effects of the fiberglass layer and the 
measurement uncertainties in the experiment as an indicator of surface azimuthal direction. 
Although the periodic water surface is an extremely simplified model of an actual ocean 
surface, this experiment further strengthens the idea of using passive polarimetry to infer 
wind direction over the ocean. Further research into this area and into other applications 
of passive polarimetric remote sensing will continue. 
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Appendix A: EBC Coupled Matrix Equation 


The matrix equation of the EBC can be written as: 
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(107) 



where ( Sj ) indicates that the integration is over surface j, 

Ooin = $no\/FoEo(ei • y) (108) 


and 

ooil = ^noy/fio~(h xii-y) (109) 

The above a 01 elements are determined by the y component of the incident electric and 
magnetic fields. Truncating the Fourier series of the surface field unknowns results in a 
finite matrix which can be inverted to obtain the surface field unknowns. 


For the special case of two surfaces described by 

fi{ x ) = h\ cos (~^~ ) (110) 

f 2 {x) = — h 2 cos (-jj-) — d (111) 

the matrix elements can be integrated and expressed in terms of Bessel functions as follows 
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[*?L “ (HD 


where «7| m _ n | denotes the Bessel function. These Bessel functions were evaluated using a 
series form, which was summed until convergence to within a fractional accuracy of 10 -6 
was observed. 
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Appendix B: MoM Matrix Equation 


The matrix equation for the Method of Moments can be written as: 
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with the testing points and integration on surface one, 
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[i] = r Gip(t m ,x)dx (130) 
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with the testing points and integration on surface two. The above integrals were approxi- 
mated by a single point rectangular area at the testing point. This approximation becomes 
more accurate as the number of testing points is increased. In the results presented, 150 
basis functions were used, which corresponds to 66 basis functions per wavelength in the 
free space region. This large number of basis functions was needed due to the much shorter 
wavelength in the high permittivity water region. 
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Table 1: Results of measurements at X band for 0 = 20' 
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Table 2: Results of measurements at X band for 6 = 30° 
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9 Figure Captions 


Figure 1: Geometry of a “two-layer” periodic surface 


Figure 2: Comparison of EBC and MOM brightness temperatures at polar angle 20° versus 
azimuthal angle (a) Tbh (b) Tjj u (c) Ub (d) Vb 


Figure 3: Comparison of theoretical and experimental brightness temperatures at polar 
angle 20° versus azimuthal angle (a) Tbh (b) Tb v (c) Ub (d) Tb p 


Figure 4: Comparison of theoretical and experimental brightness temperatures at polar 
angle 30° versus azimuthal angle (a) Tbh (b) Tb v (c) Ub (d) Tb p 
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